Astrophysical J. (in press Dec2007) 



WATER IN COMET C/2003 K4 (LINEAR) WITH SPITZER 

CHARLES E. WOODWARD 1 , MICHAEL S. KELLEY 1 ' 2 , 
DOMINIQUE BOCKELEE-MORVAN 3 , 
R. D. GEHRZ 1 

ABSTRACT 

We present sensitive 5.5-7.6 yum spectra of comet C/2003 K4 (LINEAR) ob- 
tained on 16 July 2004 (r h = 1.760 AU, A Spi tzer = 1-409 AU, phase angle 35.4°) 
with the Spitzer Space Telescope. The v-2. vibrational band of water is detected 
with a high signal-to-noise ratio ( > 50). Model fitting to the best spectrum yields 
a water ortho-to-para ratio of 2.47 ± 0.27, which corresponds to a spin temper- 
ature of 28.5115 K. Spectra acquired at different offset positions show that the 
rotational temperature decreases with increasing distance from the nucleus, which 
is consistent with evolution from thermal to fluorescence equilibrium. The in- 
ferred water production rate is (2.43 ± 0.25) x 10 29 molec. s -1 . The spectra do 
not show any evidence for emission from PAHs and carbonate minerals, in con- 
trast to results reported for comets 9P/Tempel 1 and C/1995 01 (Hale-Bopp). 
However, residual emission is observed near 7.3 /im the origin of which remains 
unidentified. 
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INTRODUCTION 



The composition of cometary nuclei probes the physical conditions in the early solar neb- 
ula, the survival of materials from the interstellar m edium (ISM), and the cold dense molec- 



ular c loud core in which t he solar system formed (IWooden et al.l 12004 ; lEhrenfreund et al. 



2004 ; iMumma et al.l 120031 ) . Comet nuclei are highly porous agglomerates of ice and dust 
grains, perha ps with highly strat i fied, inhomogeneou s layers of varied density, porosity, and 



Prialnik et al. 



composition (Harker et al.l 120071 : Belton et al.l 120061 : lOro et al.l 12009 : lA'Hearn et all 12005 



20041 ). The nucleus composition is dominated by ices (primarily water), or- 



ganic refractory materials, silicates, and carbonaceous materials. When comets are within 
heliocentric distances of < 20 AU, solar insolation triggers subl imation and the release of 
volatile gases, sometimes sporadically, forming observable comae (IMeech fc Sovernll2004l ). 



In the nucle us of a come t , volatiles are frozen as ices or trapped as gases in amor- 
phous water ice (jCaprial 120021 : iPrialnikj l2002l ). Cometary activity occurs when gases are 
released through sublimation or through the exoergic crystallization of amorphous water ice. 
Between ~ 20 to 5 AU, when nuclear surface temperature s reach ~ 20 — 100 K, CO ice 
sublimes from the nucleu s ( jCapria et al. 2000 : Prialnik 2002 ). possibly from near the surface 
( Gunnarsson et ablbood ). and triggers activity and intermittent outbursts. Between ~ 6 to 
~ 4 AU, a dramatic increase occurs in gas production and grain entrai nment and signals th e 
coma onset stage. At nuclear surface temperatures of ~ 120 — 130 K (jPrialnik et al.ll2004l ). 
the water ice phase transition (amorphous to crystalline) releases a fraction of the trapped 
volatile gases. Strong erosion maintains the CO-ice sublimation and phase transition fronts 
relatively close to the surface (jCapria et al.l 120001 ) . At ~ 4 to ~ 3 AU, the near-surface 
crystalline water ice layer, with its remaining trapped gases, begins to sublime. Water subli- 
mation drives this vigorous activity stage that is often characterized by discrete active areas 
or jets. 

Water is the dominate ice in comet nuclei and the production rate of water is correlated 
with comet activity. It influences the thermal balance of the coma as a strong coolant. At 
some wavelengths < 10 /xm emiss ion from ro-vibrationa l transitions of water can dominate 
the spectral energy distribution (ICrovisier et al.l Il997bl ): water can also be observed from 



its ro tational transitions at submillimeter wavelengths (see review of iBockelee-Morvan et al. 
20041 ). Probed through spectroscopic observations of coma species, the water production 



rate, coma temperature, and the nuclear spin temperatures derived from ortho-to-para ra- 
tios (OPR) are of particular interest in the study of cometary atmospheres and cometary 
physics. These physical characteristics, complemented by knowledge of the nucleus refrac 



tory and ice composition, p rovide constraints on solar nebula models (IMumma et al.l 12003 



Markwick fc Charnleyll2005l ). and restrict the formation zones within the protoplanetary disk 
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where cometary nuclei could conglomerate. In particular, the nuclear spin temperature of 
water measured in comet coma e may be indicative of the chemical formation temperature of 
water (jDello Russo et al.l 120051 ; iMumma et al.l Il993l ) therefore identifying the environment 
where pre-cometary ices condensed. 

Here we present longslit Spitzer Space Telescope spectroscopic observations of the 6 jum 
z/ 2 vibrational band of water detected in comet C/2003 K4 (LINEAR) at = 1.760 AU. 
The high signal-to-noise and the Infrared Spectrograph (IRS) longslit enable us to extract 
spatially resolved spectra in the coma and to measure the water production rate, Q(H 2 0), 
and the rotational temperature, T rot , and the OPR variation in the coma. Space observa- 
tions of the strong i/ 2 fundamental bands near 6 /zm present a potentially more advantageous 
method for constraining water production rates and T rot in comets than the more common 
ground-based measurement of the weaker non-resonance fluorescent "hot-bands" near 2.9 fim 
as the complex corrections for telluric extinction, slit-loss due to seeing, and con sideration of 



whether th e local radiative pump in the coma is optically th ick are minimized (jBonev et al. 



20071 . 120061 ; iDello Russo et al.ll2004l ; iBockelee-Morvanl Il987t ). In addition, accurate labora- 
tory measurements of the absorption line strengths used to compute Einstein coefficients 



for the v<i pump from the ground-state (000) are extant (IBarber et al.l 12006 



Dello Russo et al.ll2004l ; iPartridge & Schwenkdll997l ) while those for the hot bands are more 
challenging, leading to some uncertainty in estimates of the spontaneous emission rates, 

The infrared z/ 2 band of water was first detected with the Short Wavelength Spectrometer 
(SWS) of the Infrared Sp ace Observatory (ISO) in the exceptional comet C/1995 01 (Hale- 
Bopp) at r/i ~ 2.9 AU (jCrovisier et al.l Il997al ). The SWS spectral resolution of ~ 1000 
resulted in the detection of several individual ro-vibrational lines. However, the low signal- 
to-noise ratio prevented detailed analysis of their relative intensities. We discuss the Spitzer 
observations and data reduction techniques in § [23 Section I3TT1 discusses the modeling of the 
z/ 2 water band. Section 13.21 through § 13.51 present the results, followed in § 0] by a discussion 
of residual emission features, including a comparison with the z/ 2 band of water detected in 
other comets with Spitzer. Section [5] presents a summary of our study of comet C/2004 K4 
(LINEAR). 



2. OBSERVATIONS AND REDUCTION 



First identified as an asteroidal object in the LI NEAR survey, C/2 03 K4 (LINEAR) was 
discovered to have an extended spherical coma by lYoung fc McGahal (120041 ) with parabolic 
orbital elements consistent with that of a dynamically new Oort Cloud comet. The following 
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year upon perihelion approach (q = 1.02 AU, 2004 Oct 13.5 UT), C/2003 K4 (LINEAR) was 
bright in the optical (V < 7 mag) and was noted to exhibit a primarily featureless 10 jum 
spectral energy distribution with emission fro m large amorphous carb on and silicate grains 
(grain radii > 0.7 /im) dominating the coma (IWoodward et al.ll2004f). The 1 /im silicate 
featu re-to-continuum ratio was observed to be near unity (jsitko et al. 2004 ; Russell et al. 



20041 ). with lit tle evidence for structure near 11.2 jum at tributable to Mg-rich crysta lline 



olivine grains (ISitko et all 12004! ; IWoodward et al.ll2004h . ISchulz. Stiiwe. k Erdl (120051 ) ob- 
served a single broad coma feature in broadband images, perpendicular to the sunward 
direction in C/2003 K4 (LINEAR) during the 2004 May through 2004 July period with an 
increase in the amount of inferred gas-contamination in the B and V coma surface brightness 
with decreasing heliocentric distance. During the epoch of our Spitze r observations (§ 12.11) . 
the c omet produced a considerable amount of dust (Afp ~ 10, 000 cn uSchulz. Stiiwe. &; Erd 



2005|). 



2.1. Spitzer IRS 

Spectra of co met C/2003 K4 (LIN EAR) were obtained with the Infr ared Spectrograph 



(IRS) instrument (IHouck et al.l 120041 ) on the Spitzer Space Telescope ( IGehrz et al.l 120071 ; 
Werner et al.l 120041 ) . The comet was observed in the second order of the short- wavelength, 
low-resolution module (SL2) on 2004 Jul 16 at 04:56 UT as part of a Spitzer Guaranteed 
Time Observation (GTO) program (PI: R.D. Gehrz), program identification (PID) 131, 
astronomical observation request (AOR) key 0008525056, and processed with IRS reduction 
pipeline S15.3.0. The SL2 slit is 3.7" wide and provides 57" of spatially resolved spectra 
(1.8" pixel -1 ) with a spectral dispersion of 0.06 /im. Six spectra (14 s x 3 cycles) at 5.2- 
7.6 /im were recorded in a 3 x 2 spectral map, with 7.2" x 78" steps (perpendicular x parallel 
to the long slit dimension). The comet was at a heliocentric distance (r^) of 1.760 AU, a 
Spitzer-comet distance of 1.409 AU, and a phase angle of 35.4°. 

At the time of acquisition, we attempted to acquire the comet nucleus with the Spitzer 
IRS 15 /im peak-up array. However, the bright inner coma saturated a 76" x 64" ellipse in 
the 98" x 72" peak-up array preventing the spacecraft from computing a centroid on the 
comet. Thus, the telescope pointed to the comet's nominal ephemeris position derived from 
orbital elements uploaded to the spacecraft prior to the execution of the AOR. On the date of 
observation, 2004 Jul 16 UT, the nominal position of the comet derived from these elements 
was 28" from the actual position calculated with revised elements from JPL ephemeris #96 
(computed 2006 Dec 14 with a data-arc spanning from 2003 May 28 through 2006 Nov 17). 
However, comet C/2003 K4 (LINEAR) had an extensive coma (> V in diameter) at the 
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epoch of our Spitzer observation and thus error in the position of the nucleus did not affect 
our ability to obtain spectra of the comet coma. Our discussion of the IRS slit positions 
within the coma of the comet are referenced to the actual position as computed from the 
most recent JPL ephemeris. Figure [T] shows the blue peak-up image (saturated core), the 
slit positions, and the position of the nucleus (cross). 

Coma emission (including the spectral signature from water lines) is present in all por- 
tions of our slits to varying degrees. Therefore, a robust estimate of the background emission 
is difficult to accurately assess (to the level of a few percent) using the Spitzer longslit ob- 
servations of comet C/2003 K4 (LINEAR) alone. However, in the same IRS campaign 
(#10) an observation with the same IRS AOR parameters toward a similar ecliptic latitude 
(52.5° ± 0.5°) was available. The background derived from this IRS observation (AOR key 
0004733952 obtained from the Spitzer archive) was two-dimensionally subtracted from the 
basic calibrated data products (BCDs) of comet C/2003 K4 (LINEAR). 

After background subtraction, we corrected the world coordinate system of the two- 
dimensional spectral frames for the motion of the comet, then combined each s ource obser- 



vatio n into data cubes with the CUbe Builder for IRS Spectra Maps (CUBISM, ISmith et al. 
2007 ) program (vl.5)0. CUBISM combines each cycle and each slit position into a data cube 
where two axes contain the spatial information (1.85" pixel -1 grid), and the remaining axis 
contains the spectral information. A separate cube is created for the pipeline errors derived 
from the individual BCDs. The program photometrically calibrates the data, including a 
correction for diff r action losses at the entrance slit (the so-called slit-loss correction factor; 



Kelley et al.l 12006c ISpitzer Science Centerl 120061 ) . 



IRS spectra extracted from SL2 module have weak fringing artifacts ( < 2% of the 
source flux, -F(Aj), at a given wavelength) that are difficult to completely remove using 
a sinusoidal function, as they are not spectrally resolved and vary with position in the slit 
( ISpitzer Science Centeilbood ). Thus, to account for any potential residual fringe contamina- 



tion we compute the flux uncertainty in our extracted spectra at a given wavelength from the 
quadrature sum of the photometric error arising from the pipeline processing of individual 
BCDs plus a contribution due to SL2 fringing signal equal to 0.02 x -F(Aj). This latter term 
is an upper limit to the fringe uncertainty. 

We extracted spectra from nine locations in the coma, as shown in Fig. [IJ The extraction 
apertures are 1.85" x 7.40" rectangles (subtending 1890 km x 7560 km within the coma). 
We restricted our nine source extractions and subsequent analysis to locations where the 
coma is brightest, from 0" to +30" offset from the nucleus (Fig. [T]). 



1 Available at |http: / /ssc.spitzer.caltech.edu/ a rchanaly / contributed / cubism / 
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2.2. Wavelength Calibration 

The water lines are unresolved and blended in the IRS spectra. Furthermore, the Spitzer 
IRS spectra are calibrated with standards measured near the center of the slit and three of our 
extractions occur near the slit edge. Analysis of the wavelength calibration and unresolved 
line widths is critical for identifying and fitting the water lines in the IRS spectra. To 
verify the wavelength calibration at the slit edge, we reduced IRS calibration observations 
of NGC 7027 at the center and edge of the SL2 slit. The NGC 7027 observations (AOR 
key 0010066432, Spitzer PID 1410, IRS pipeline S13.2.0) were taken during the same IRS 
observing campaign (IRS #10) as the C/2003 K4 ( LINEAR) spectra. We fit Gaussians to 



the [Mg V] 5.61 /zm line ([Bernard Salas et al.l 120011 ) with the nebula at the center and edge 



positions. The width of the [Mg V] line in NGC 7027 at the center position agrees with 
the spectral resolution solution provided by the Spitzer Science Center (0.0605 jum). At the 
edge position, the line width increases to 0.0655 ±0.0002 /im. The [Mg V] line was observed 
at a central wavelength of 5.6242 ± 0.0003 /im. Accounting for Doppler shift, the observed 
wavelength is 0.0137 ± 0.0007 fim from the vacuum rest wavelength of 5.6099 ± 0.0006 /im. 
The IRS SL2 wavelength calibration is ±0.006 /im (r.m.s.), indicating the [Mg V] shift is 
2.3-cr from the expected central wavelength. The Gaussian fits show no significant difference 
in central wavelength between the center and edge positions. We compare the wavelength 
positions of the water lines to the expected central wavelengths in § 



3. THE v 2 WATER BAND IN COMET C/2003 K4 (LINEAR) 

3.1. Model Fitting 

At the resolution of the SL2 Spitzer spectrometer [R » 100), the v 2 water band shows 
ro- vibrational structure from which information on the rotational temperature, T rot , in the 
ground vibrational state can be obtained. Though the spectral resolution is not high enough 
to separate individual ro-vibrational lines, and therefore ortho from para water lines, it is 
still possible to assess whether our Spitzer spectra can provide some constraints on the 
OPR. Previous determinations of the OPR in cometary water were ba sed on water infrared 
spectra obtained with resolving powers between 1500 and 25000 (e. g., 



Crovisier et al. 1997b; Dello Russo et al. 2005; Kawakita et al. 2006; Bonev et al 



Mumma et al 



1993 



200 



7). 



Vibrational emission from cometary parent molecules results from radiative excitation 
by solar infrared radiation followed by fluorescence. For the fundamental vibrational bands 
of water, including v%, emission is not pure resonant fluorescence, as these bands are signifi- 
cantly populated by radiative decay from higher excited vibrational states. The vibrational 
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fluorescence scheme of cometary water is presented by iBockelee-Morvan fc Crovisierl (119891 ). 
The z/ 2 band is significantly populated by decay of the v 2 + v% band. The resulting emission 
rate of z/ 2 is 2.41 x 10~ 4 s _1 at = 1 AU from the Sun. The z/ 2 + z/ 3 — u 3 hot band has an 
emission rate of 7.4 x 10~ 6 s _1 and therefore does not contribute significantly to the emission 
observed between 6 and 7 //m. Other hot bands (e.g., {v\ + v 2 + z/ 3 ) — (i>\ + z/ 3 )) are even 
weaker. 



We used the model of fluorescence water emission presented by lBockelee-Morvan fc Crovisier 



(119891 ) for analyzing the Spitzer data. This model considers five excited vibrational states 



an d their subsequen t radia tive cascades, and is an improvement of that presented in detail 



in 



Bo ckelee-Morvanl (119871 ). where only the z/ 2 and z/ 3 bands are considered. Einstein co- 



efficients for ro- vibratio nal transitions are computed using the 2003 edition of the GEISA 



spectroscopic database (jJacquinet-Husson et al.ll2005l ). which includes all significant routes 



leading to v 2 excitation, including via hot bands. Comparison of the line strengths g i ven in 



GEISA with those resulting from the ab-initio calculations of iPartridge fc Schwenkd (119971 ) 
verified that the line-by-line relative intensities are insensitive (within 3-4%) to the choice 
of water line lists. 

Our water model takes into account opacity effects in vibrational excitation and emis- 
sion, using the escape probability formalism. For computing the line-by-line fluorescence, we 
considered 32 ortho and 32 para rotational levels in each vibrational state. The rotational 
populations in the ground vibrational state can be described by a Boltzmann distribution at 
a temperature T rot , or the populations can be computed using an excitation model that con- 
siders the evolving excitation conditions experienced by the water molecules as they expand 
in the coma (§ 13. 3p . In this detailed model, ro- vibrational line intensities are computed for 
a circular aperture centered on the nucleus. We do not expect the results to be significantly 
sensitive to the shape of the aperture, providing the aperture area is conserved. From the 
model output, synthetic Spitzer spectra were generated by convolving the intensity of the 
individual ro-vibrational lines with the instrumental spectral response of the spectrometer, 
described by a Gaussian. Figure [2] shows examples of synthetic spectra obtained for the spec- 
tral resolution of SL2 and a 15 times higher resolving power. At first glance, the position 
and relative intensities of the peaks in the water modeled spectrum match approximately 
those in the C/2003 K4 (LINEAR) Spitzer spectra shown in Fig. [31 thereby demonstrating 
that these spectra are dominated by water emission. 

For fitting the observed spectra, we assumed that the rotational populations of the 
ground vibrational state follow a Boltzmann distribution. The gas expansion velocity, v exp , 
was fixed to 0.8 km s . The water photodissociation rate was taken equal to 1.6 x 10~ 5 s _1 
(rh = 1 AU). This rate takes into account the solar activity at the time of the observations 



following the formalism described by lCrovisierl (119891 ) . The only free parameters of the model 
are the water production rate Q(H 2 0), T rot , and the OPR. 



Overall, opacity effects are small. We computed that they affect the total intensity of 
the u 2 band by 6% for the spectrum of C/2003 K4 (LINEAR) acquired closest to the nucleus 
(7.2" offset, Fig. [p. In addition, if opacity effects are not properly taken into account in the 
calculations, then the derived OPR also can be underestimated (on the order of ~ 6%). 

The water band emission is in excess of the dust continuum emission (Fig. [3]). In the 
first analyses, the underlying continuum was determined using a 5-th order polynomial fit, 
and the residual (continuum-subtracted) spectra were fit with the water model, applying 
a least-squares method that uses the gradient-search algorithm of Marquardt. Continuum 
subtraction was not completely satisfactory, as excess continuum emission remained near 
6.26 /im, while the u 2 band is almost free of lines at this wavelength (Fig. [2]). More robust 
fits could be obtained by fitting simultaneously the underlying continuum and the water 
emission. Thus, we fit the original spectra with a composite curve consisting of the mod- 
eled water spectrum superimposed on a polynomial. Polynomials of 5— or 6— degree were 
used. However, final results were not found to be significantly sensitive to the choice of the 
polynomial degree between 3 and 6. 

The nominal spectral resolution of SL2 is 0.060 fim near 6 /im. Model fits to the 
best C/2003 K4 (LINEAR) spectrum (7.2" offset) with the spectral resolution left as a free 
parameter yielded AA = 0.067 ±0.004 /im, agreeing with the edge observation of NGC 7027 
(§ 12.21) . Results of model fits given in Table [1] were obtained with AA fixed to 0.065 /im. 
However, almost identical results (within the error) are returned with AA = 0.060 /im. For 
example, for the highest signal-to- noise ratio spectrum (Fig. [3l box A), the retrieved OPR 
is changed from 2.47 ± 0.27 to 2.31 ± 0.24. We also found that the frequency calibration in 
SL2 spectra is likely incorrect by a tenth of the resolving power. The central wavelengths 
of v 2 band structures in C/2003 K4 (LINEAR) spectra are better matched by shifting the 
observed spectra by 0.0032 to 0.0062 /im, within the errors discussed in § [2j For example, 
for the 7.2" offset spectrum, the \ 2 between 5.8 and 7.1 jum is decreased by a factor 2.8 when 
applying a 0.0062 /im offset. The spectra of comet C/2003 K4 were shifted by 0.0032 to 
0.0062 /im for the model fits shown in Table [Hand the corresponding figures (Figs. [3] and H]). 



3.2. Water Modeling Results 



The best-fit modeled spectra for comet C/2003 K4 (LINEAR) are shown in Figs. [3] 
and HI where in the latter figure the continuum has been subtracted. Residuals with respect 
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to observed spectra are shown in the bottom of Fig. HI Retrieved model parameters are given 
in Table [TJ The agreement between our models and the Spitzer spectra is rather good for 
slit extractions A, D, and G, with reduced y 2 less than 1 (Tabled]). 

Some excess emission (> 3-cr deviation) is noticeable at 6.05 /im in most of the spec- 
tra at offset > 18". Models that incorporate values of T rot higher than those determined 
from model fitting (Table [1]) reduce the residual continuum emission near 6.05 //m but are 
inconsistent with the relative water line intensities arising from intrinsically stronger lines 
(Fig. [2]) measured longward 6.3 /zm. The variation in emergent water line intensities for 
20 < T rot (K) < 90 which produce the broad emission feature from 5.5 to 7.0 /im when 
the models are convolved to the resolution of the Spitzer IRS SL2 are shown in Fig. [51 The 
6.05 /mi peak (mainly ortho 2 12 -l i line) is more intense than the 6.18 /im (mainly l 10 -loi 
line) peak only for high T rot (> 60 K). At low T rot , these two lines result essentially from IR 
pumping from the loi ground state rotational level: the ratio of their intensities 7(6.18 /xm)/ 
7(6.05 /im) then depends uniquely on ro- vibrational Einstein ^-coefficients and is predicted 
to be 1.5. 

Weak residual emission (2-er deviation) is also observed at 5.90-5.95 /im (see B, E, H 
extractions in Fig. [3]). This excess emission does not coincide in wavelength to the position of 
the 5.88 /xm peak of the water band, and may be related to flaws in background subtraction. 

Our derived values of T rot and OPR could be affected by the residual emission present 
between 5.90-6.1 /im. This emission (in excess of fringe artifact contamination) may arise 
from sources of weak line and continuum emission other than water not accounted for in our 
models. To quantify such effects, we derived estimates of T rot and OPR by independently 
fitting two partial subsections of the spectra (5.85-6.3 jum and 6.3-7.0 /im) and any con- 
tinua shortward of 5.8 /im and longward of 7.0 /im for the brightest slit extractions. These 
independent model fits, including those derived from the best reduced x-squared fit of the 
entire spectra, are summarized in Table [2j Consistent OPR values are obtained from the 
independent model fits. The derived T rot values are higher as a result of the 5.90-6.1 /im 
excess flux only when the 5.85-6.3 /im part of the spectrum is considered. Table [1] provides 
model parameters for all slit extractions retrieved by fitting the 6.3-7.0 /im part of the water 
spectrum. However, these latter fit parameters are not significantly different from those 
derived by fitting the entire spectrum. 
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3.3. Rotational Temperature 



The inferred rotational temperatures are between 18 and 34 K and tend to decrease 
with increasing distance from the nucleus (Fig. [6]). The rotational temperatures are not 
constant with respect to cometocentric distance out to 3 x 10 4 km within the uncertainties. 
A linear least squares fit to T rot in Table HI as a function of distance to nucleus, gives a 
slope of —(6.1 ± 2.2) x 10~ 4 K km -1 . Such behavior is expected. In the inner coma, 
collisions are important and the rotational levels of the fundamental vibrational state are 
thermalized at the local temperature. However, in the outer coma radiative pumping prevails 
and t he rotational population of the water molecules reaches a cold fluorescence equilibrium 



^e.g., iBockelee-Morvanl 119871 ). How the populations evolve from thermal to fluorescence 
equilibrium depends of the density of the collisional partners, and related collisional cross- 
sections. Other attemp ts to examine the variation in water rotational temperature with 
cometocentric distance (jBonev et al. 2007 ) are restricted to the inner ( < 10 3 km) collision- 
dominated coma. 

Water excitation model s currently developed by various investigators include both H 2 0- 
H 2 and H 2 0-e" collisions jBiverl Il997l : iBensch k Berginlhooi IZakharov et al.ll2007h . For 



our analysis, we use the exci tation model of 



Bensch fc Berginl ( 120041 ) and IZakharov et al. 



Biverl (119971 ). which differs from the models of 



(120071) by the method use d to solve radiation 



trapping effects, but yields almost similar results (IZakharov et al.ll2007l ). The electron den- 
sity and radial temperature distribution is based on the measurements o f lP/Halley made 



by the Giotto mass spectrometers (e.g., lEberhardt fc Krankowskyl Il995l ). to which scaling 



factors are applied to account for variations with water production rate and heliocentric dis- 
tance. The parameter x ne is a multiplying factor to the electron density, normalized to the 
lP/Halley Giotto measurements (x ne (lP/Halley) = 1). The output of the model is the rota- 
tional populations in the ground vibrational state as a function of distance to nucleus. The 
populations are included in the water infrared fluorescence model to simulate water spectra 
at offset positions. For direct comparison with the observations, rotational temperatures are 
derived by fitting the synthetic spectra, as was done for the observed spectra. 

Figure [6] shows the evolution of T rot predicted for x ne values of 0.2, 0.5, and 1.0, and inner 
coma kinetic temperatures Tki n of 30, 40, 50, and 100 K. These values of the x ne parameter 
were selected because the el ectron density is r ather uncertain and mapping of the 557 GHz 
H 2 line favors x ne ~ 0.2 (IBiver et all 120071 ). The predicted increase in T rot at ~ 2000 
km offset is due to thermal excitation by hot electrons. In the electron density model, this 
distance corresponds to t he contact surface, Re s, where the electron temperature and density 
undergo a steep increase ( Xie fc Mumma 19921 ). Beyond ~ 2Rcs, T rot decreases because the 
effect of excitation by electronic collisions become less efficient with respect to radiative decay 
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and the fluorescence equilibrium of the ground vibrational state is cold. In the 7000-15000 
km region, the model predicts a decrease in T rot of AT rot = 6 to 11 K for the considered 
parameters, in contrast to the observed decrease of ~ 5 to 6 K. If one considers uncertainties 
associated with T rot derived from the Spitzer spectra, models with T kin = 30 — 90 K and 
x ne = 0.2 — 0.5 are satisfactory. A x ne value of 1.0 d oes not fit the data obtained at 7.2" 
offset (Fig. [6]), in agreement with iBiver et al.l (120071 ) . The kinetic temperature is poorly 
constrained at the sampled cometocentric distances because, T rot retains little memory of 
excitation conditions prevailing in the collisional region. 



3.4. Ortho-to-Para Ratio 

Though the spectral resolution of the SL2 spectra is low (R ~ 100), model fitting pro- 
vides an accurate measurement of the OPR in comet C/2003 K4 (LINEAR) for the high 
signal-to-noise spectra (Tables [T] and [2]). The OPR can be retrieved from u 2 band spectra 
obtained at low resolving power because several para lines are well separated in wavelength 
from strong ortho lines (Fig. [2]). The band regions most sensitive to the OPR lie at 6.12 
and 6.4 /im. At 6.12 /im, emission is dominated by the ln-Ooo para line, and three other 
significant para lines (Fig. Since the nearby 6.05 and 6.18 /xm peaks are mainly due 
to ortho lines (2 12 -l i and l 10 -l i, respectively), the intensity ratios 7(6.12/xm)//(6.18/xm) 
and J(6.12/im)//(6.05/xm) increase with decreasing OPR. Similarly, at 6.4 /xm, the contri- 
butions from para lines (Ooo~lii and ln-2oa) dominate the spectrum and the intensity ratio 
I(6.4/xm)//(6.5/xm) is a function of the OPR. The variation in water v-i band features with 
OPR, for T rot = 30 K, at the Spitzer IRS spectral resolution is illustrated in Fig. From 
the models depicted in this figure, the peak intensities of the 6.05 /xm, 6.18 /xm, and the 
6.50 /xm features increase by ~ 8%, ~ 8%, and ~ 13% respectively as the OPR changes from 
values of 2 to 3. The intensity at 6.12 and 6.4 /xm decreases in turn by ~ 18%. However, 
the 6.64 and 6.85 /xm water features remain constant and their intensity ratio only depends 
upon T rot . The different behaviors of the intensity of the features with T rot and OPR make 
the accurate measurement of these two parameters possible. 

For the comet C/2003 K4 (LINEAR) spectrum at position A, 7.2" offset from the 
nucleus position (Fig. [1]), we obtain an OPR = 2.47 ± 0.27 when fitting the entire water 
spectrum (Table [2]). The other spectra yield OPRs consistent with this value (Tables [Tj, Ui>- 
Conversions between ortho to para states by radiative transitions or by collisions in the coma 
have very low probability. The constancy of the O PR in the coma also has been convincingly 



demonstrated in comet C/2004 Q2 (Machholz) by iBonev et al.l (120071 ). The weighted mean 



of all OPR values in comet C/2003 K4 (LINEAR) given in Table H (fits to 5.8-7.0 /im region) 
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is 2.43 ± 0.15. The OPR value derived for the aperture slit closest to the nucleus (labeled 
A), 2.47 ± 0.27, corresponds to a spin temperature T spin = 28. K. 



'-3.5 



The reduced x 2 between 5.8 and 7.1 /mi obtained for the spectrum at 7.2" offset is 0.5. 
When the OPR is fixed to OPR = 3, the reduced x 2 is — 15% higher (45% higher when 
fitting 6.3-7 fim partial spectrum). Figure [8] shows the model fit obtained in this case which 
yields T rot = 30.5 ± 3.3 K, a value close to that obtained with OPR = 2.47 (Table [2]). There 
is significantly higher discrepancy between 6.3 and 6.4 /im in the two models. 

The OPR (and T rot ) determination relies on the assumption that water emission domi- 
nates the 5.8-7.2 /mi C/2003 K4 (LINEAR) spectrum. Misleading results can be obtained 
when extra emission from other constituents is prese nt. In addition to the 6 /zm residual 
emission discussed in § 13.21 PAH emission peaks (e.g.. iPeeters et al.ll2002l ) near 6.2 /im, i.e., 
in the region where the shape of H 2 band depends on the OPR (Fig. [8]). Because the 
independent model fits of partial subsections of the spectra provide consistent OPR val- 
ues (Table [2]), our OPR determinations are likely not affected by unaccounted species. We 
also conclude that the dominant source of emission in the SL2 wavelength regime for comet 
C/2003 K4 (LINEAR) is the water u 2 band. 

Our derived water OPR for C/2003 K4 (LINEAR) is comparable to to values derived 
for Oort Cloud (nearly-isotropic) comets such as C/1995 01 (Hale-Bopp), C/1999 HI (Lee), 
C/1999 S4 (LINEAR), C/2001 Ql (NEAT), and C/2001 A2 (LINEAR), although lower (~ 
15-20%) than that reported for Jupiter-family (ecliptic) comets, for example 103P /Hartley 2 



or lP/Halley (ICrovisier et al. 



Q2 (Machholz) flBonev et al 



2000; lDello Russo et al.ll2005l : iKawakita et alJl2006h or C/2004 



20071 ). The value of T t 



spin 



28.5Z|:| K for C/2003 K4 (LIN- 
EAR) is suggestive of precometary ice formation in a cold molecular clou d environment 
devoid of secondary processing in a warm solar nebula (IKawakita et al.ll2006l ). although the 
precise interpretation of the OPR as a probe of the pr imordial formation zones of comets in 
the protosolar nebula remains vexing (jCrovisierl 120071 ) . We also do not have spectroscopic 
measurements of other common cometary ices such as ammonia or methane in coma of 
C/2003 K4 (LINEAR). Thus we are unable to ascertain whether the ices incorporated into 
the nucleus of C/2003 K4 (LINEAR) share the same chemical compo sition and hom ogeneity 



of T spin as that found for other comets of diverse dynamical classes (jCrovisierl 120071 ) 



3.5. Water Production Rate 



The intensity of the v-i band measured for C/2003 K4 along the nine slit extractions 
is given in Table [TJ, and is plotted as a function of offset in Fig. The evolution with 
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distance to nucleus is consistent with that computed using a Haser distribution for the water 
density and Q(H 2 0) = (2.43 ± 0.25) x 10 29 molec. s" 1 , where th e error includes a 10% 
uncertainty in the IRS absolute calibration (jSpitzer Science Center! 120061 ). Some deviations 
are observed, which may be related to asymmetries in the density distribution and/or (for 
the noisy spectra) incorrect background subtraction. 

The derived pre-perihelion (r^ = 1.76 AU) water production rate is consistent with OH 
18-cm observations performed with the Nangay radio telescope which yield Q(H 2 0) ~ 2 x 
10 29 molec. s _1 at the epoch of the Spitzer observations (Crovisier et al., personal communi- 
cation). Post-perihelion measurements obtained from H2O 557 GHz line observ ations using 
the O din satellite give a Q(H 2 0) about a factor of two lower at = 1.7—1.8 AU (IBiver et al. 



20071 ). which suggests a pre/post-perihelion asymmetry in the gaseous activity of the comet. 
Similar asymmetrical perihelion production rates of water and other vola tiles has been ob - 
served in other comets, including C/1 995 01 (Hale-Bopp ) (e.g., Fig. 3 of IBiver et al.lll997l ) 
and possibly lP/Halley (e.g., Fig. 6 of iGehrz et al.l 1200 



4. Residual Emission Between 5.5 and 7.6 /im 



Lisse et al.l (120061 ) report emission from carbonate minerals at 6.5-7.2 /im and organic 



(PAH) emission at 6.2 jum in the spectrum of comet 9 P /T e mp el after collision with the 
Deep Impact impactor. In a subsequent paper, iLisse et al.l (120071 ) claim the detection of 
these emission f eatures in the ISO spectrum of comet C/1995 01 (Hale-Bopp) published by 
Crovisier et al.l (1997b) . However , a re-a nalysis of the ISO observations of comet Hale-Bopp 



Crovisier fc Bockelee-Morvanl (120071) does not confirm the det e ction of PAHs reported 



Lisse et al.l (120071 ). Furthermore, ICrovisier fc Bockelee-Morvanl (120071 ) demonstrate that 



by 
by 

carbonate emission at 7 /xm, though possibly marginally present, is fainter by a factor 2 to 
3 than asserted by ILisse et al.l (120071 ). Figure [8] shows representative spectra of PAHs and 
carbonate minerals compared to the best Spitzer spectrum of comet C/2004 K4 (LINEAR). 
Our synthetic water spectrum wholly accounts for any emission in-excess of the continuum 
at the wavelengths of PAHs and carbonate emission. As discussed in § 13. H and § 13.21 residual 
emission is only marginally present in some spectra near 5.9 and 6.05 /im in various aper- 
tures in the C/2003 K4 (LINEAR). Since PAH features are narrow (AA ~ 0.15 /im) and 
peak near the water 6.18 /im pattern, a significant contribution of PAHs in the spectrum 
would have resulted in an intensity ratio 1(6.18 /iin)//(6.05 /im) higher than observed. Sim- 
ilarly, carbonate emission, if present, would have been seen directly on the original spectra 
(continuum background included, Fig. [3]) longward of 7 /mi. Indeed the 3-<r upper limit to 
the peak intensity of any carbonate or PAH emission, computed from the residual emission 
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between 6.15 to 6.30 jum and 6.75 to 7.25 /im using the representative PAH and carbonate 
spectra shown in Fig. [8], does not exceed 7 to 8 x 10 -21 W cm -2 /im" _1 ( < 10 -3 Jy) in the 
C/2003 K4 (LINEAR) spectrum at 7.2" offset position. 

As depicted in Fig. M, Spitzer spectra of comets C/2004 Bl (LINEAR), 71P/Clark, 
and 9P/Tempel 1 show evidence for i/ 2 water emission arising from sublimating ices in their 
comae. A detailed study, similar to that presented in this paper, i s required to investigat e 
whether emission from other compounds is present in these spectra (IWoodward et al.ll2007l ). 



From Fig. [HI we see that a narrow (AA ~ 0.1 /xm) residual emission feature is present 
near ~7.3 /im (1370 cm -1 ). A small spectral segment of the Spitzer IRS data near 7.3 /im, 
shown in Fig. [TTJ, provides detailed, close-up view of this emission feature. The origin of this 
feature is unclear. The peak flux and integrated feature flux is in excess of that anticipated 
from spurious fringe signal power. A possible candidate is emission from the SO2 ^3 band 
at 7.34 ixm that has a fluorescence emission rate at 1 AU from the Sun of 6.6 x 10 -4 s _1 
(ICrovisierl I2002T) . Syn thetic spectra of SO2 ^3 band obtained using the HITRAN database 



(IRothman et al.l 120041 ) approximately match the width of the feature, but the central wave- 



lengths do not coincide. In addition, the measured intensity in the spectrum obtained at 7.2" 
offset in the coma of comet C/2003 K4 (LINEAR), ~ 10~ 21 W cm" 2 , would imply a S0 2 /H 2 
production ra te ratio of 2.5%, a factor 10 t imes higher than measured in comet C/1995 01 
(Hale-Bopp) (IBockelee-Morvan et al.ll2000l ). Therefore, it seems unlikely that the observed 
7.3 yum feature is due to S0 2 . The NIST databas^l provides band positions of a number of 
gas phase species, including organics. No satisfactory candidate cou ld be found. For exam- 
ple, m ethyl formate HCOOCH3, identified in cometary atmospheres (IBockelee-Morvan et al. 
20001 ). has a band of medium strength at 1371 cm -1 , but also a much stronger band at 1754 



cm -1 (5.7 /im) which is not seen in the Spitzer spectra. Acetic acid exhibits a strong band 
at 1375 cm -1 , but a still stronger one at 1248 cm -1 (6.94 fim). If the feature is originating 
from a gas phase species, then the abundance of this molecule relative to water should be on 
the order of 1% or more, based on the measured intensity and typical fluorescence emission 
rates in cometary environment. 

The 7.3 fim wavelength corresponds to the characteristic vibrational frequency of the 
CH 3 "umbrella" deformation mode (~1375 cm -1 ). A 7.3 1 1m absorption feature has been 



detec ted in some galactic and extragalactic sources (e.g., IChiar et al.l |2000| ; ISpoon et al. 



20001 ). and assigned to aliphatic hydrocarbons. However, an absorption signature at 6.85 /xm 
is also observed, corresponding to CH 2 bending vibrations. Spectra of various carbonaceous 
refractory materials, including chondritic material, that contain aliphatic chains show that 



2 http: / /webbook. nist.gov/chemistry/ 
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these two features are present a s a doublet with intensity ratio /(6.8/mi)//(7.3/im) > 1 
( IPendleton &: Allamandolall2002l ). In contrast, no residual emission is observed at 6.8 /im in 
the comet C/2003 K4 (LINEAR) spectrum. Therefore, aliphatic hydrocarbons are not likely 
the source of the cometary 7.3 /im feature, although such compounds have be en identified 



in th e material captured from comet 81P/Wild 2 by the Stardust spacecraft (jKeller et al. 
2006h . 



A few galactic and extragalactic sources exhibit a weak emission feature at 7.3-7.4 /im (with 
no 6.85 /mi counterpart) that shows up on the wing of the well known strong 7.7 /mi complex 
attribut ed to CC stretching/ CH in plane bending vibrations of aromatics (likely PAHs) com- 
pounds (IPeeters et al.ll2002l ). Based on theoretical calculations of exp ected CC stretching 
band positions of PAHs of various complexity (see IPeeters et al.ll2002l ). this weak compo- 
nent is likely a PAH signature. However, the comet feature may have a different origin as 
no strong 6.8 and 7.7 /im PAH emission is evident in the Spitzer spectra of C/2003 K4 
(LINEAR). 



5. CONCLUSIONS 

We have observed the z/ 2 vibrational band of water in comet C/2003 K4 (LINEAR) 
within 5.5 to 7.6 /tm Spitzer IRS spectra, deriving a water production rate of (2.43 ±0.25) x 
10 29 molec. s -1 when the comet was at a pre-perihelion heliocentric distance of 1.760 AU. 
Although the IRS spectra are of moderate resolution, modeling of the observed emission 
in the 5.7 to 6.8 /tm region constrained the water spin temperature to be 28.5l 3 J K. The 
measured T spin is comparable to that of other Oort cloud comets and suggestive of a common 
formation zone for the precometary water ices that eventually agglomerated into the nuclei, 
though the precise interpretation of the OPR as a probe of the primordial formation zones 



of comets in the protosolar nebula remains controversial (ICrovisierl 120071 ) . The observed 
decrease (at 3-cr confidence level) of the water rotational temperature with cometocentric 
distance is compatible with evolution from thermal to fluorescence equilibrium and constrains 
somewhat the role of electron collisions in water excitation. The kinetic temperature of the 
gas is poorly constrained. 

Neither emission from carbonates nor PAHs was necessary to account for any emission 
in excess of the continuum at wavelengths between 5 to 7 /tm, suggesting that these species 
are not presen t in t he coma of C/2003 K4 (LINEAR) at the abundance levels measured 



by iLisse et al.l (120061 ) in comet 9P/Tempel 1. However, an emission feature at ~ 7.3 /im is 
observed that remains unidentified, as potential emission candidates, the SO2 ^3 band or 
CH 3 deformation modes, can be discounted. 
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Table 1. WATER IN C/2003 K4 (LINEAR): MODEL FITS 



Extraction 


Offset* 


Offset a 


Band Intensity 13 


J-rot 




OPR c 


~.2c,d 


Slit ID 


(") 


(km) 


(10~ 20 W cm" 2 ) 


{K) 








A 0-27 


7.2 


7377 


5.51 ± 0.16 


33.5 ± 


4.0 


2.47 ± 0.30 


0.5 


B 0-22 


12.9 


13198 


2.67 ± 0.10 


28.6 ± 


3.6 


2.47 ± 0.34 


1.0 


C 0-17 


21.3 


21785 


1.62 ± 0.09 


22.8 ± 


4.3 


2.5 C 


1.6 


D 1-27 


14.4 


14682 


2.88 ± 0.10 


27.6 ± 


3.5 


2.13 ± 0.31 


0.5 


E 1-22 


18.0 


18348 


1.93 ± 0.09 


20.9 ± 


3.5 


2.56 ± 0.47 


1.5 


F 1-17 


24.7 


25229 


1.19 ± 0.09 


24.8 ± 


5.6 


2.5 C 


1.6 


G 2-27 


21.7 


22154 


1.68 ± 0.08 


22.4 ± 


4.2 


2.61 ± 0.31 


0.9 


H 2-22 


24.2 


24750 


1.46 ± 0.08 


23.2 ± 


4.3 


2.38 ± 0.48 


1.5 


I 2-17 


29.6 


30209 


1.02 ± 0.09 


17.7 ± 


5.0 


2.5 C 


2.4 



a Relative to the comet C/2003 K4 (LINEAR) ephemeris position derived from 
JPL #96. 

b Band intensity above fitted continuum between 5.8 and 7.1 /xm. 
c Fits to the 6.3-7.0 /im region. 

d Reduced x 2 between 6.3 and 7.0 fxm (xl — X 2 /dof, with number of degrees of 
freedom [dof] = 23-3 = 20). 

c Assumed value. 
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Table 2. WATER IN C/2003 K4 (LINEAR): MODEL FITS USING SEGMENTS 



Extraction 


Regions Fit 


T 

-*■ rot 


OPR 


Slit ID 


(/im) 


(K) 




A 0-27 


5.85 


- 7.00 


30.7 ± 3.2 


2.47 ± 0.27 


A 0-27 


5.85 


- 6.30 


36.5 ± 5.8 


2.03 ± 0.40 


A 0-27 


6.30 


- 7.00 


33.5 ± 4.0 


2.47 ± 0.30 


B 0-22 


5.85 


- 7.00 


26.5 ± 3.3 


2.95 ± 0.40 


B 0-22 


5.85 


- 6.30 


59.5 ± 8.8 


2.83 ± 0.68 


B 0-22 


6.30 


- 7.00 


28.6 ± 3.6 


2.47 ± 0.34 


D 1-27 


5.85 


- 7.00 


27.0 ± 2.9 


2.07 ± 0.24 


D 1-27 


5.85 


- 6.30 


42.7 ± 7.6 


1.55 ± 0.30 


D 1-27 


6.30 


- 7.00 


27.6 ± 3.5 


2.13 ± 0.31 



G 2-27 
G 2-27 
G 2-27 



5.85 - 7.00 21.4 ± 3.2 2.67 ± 0.41 
5.85 - 6.30 54 ±12 2.16 ± 0.57 
6.30 - 7.00 22.4 ± 4.2 2.61 ± 0.31 
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Offset (arcsec) 



Fig. 1. — Left: The 15 /im IRS acquisition image and observed slit positions in comet 
C/2003 K4 (LINEAR). The dotted contour outlines the saturation ellipse in the image and 
the dashed line traces the edge of the peak-up array. The JPL #96 ephemeris position of 
the comet nucleus is marked with a cross (see § [2]). The slit positions are outlined with 
solid lines. Celestial north (N) and east (E), the comet's heliocentric velocity (v), and the 
direction of the sun (©) are marked with arrows. Right: The nine slit extraction apertures 
are outlined with solid lines and run from left to right, and top to bottom (see Tabled]) A, 
B, C (row 1), D, E, F (row 2), and G, H, I (row 3). 
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Fig. 2. — Synthetic spectra of the water z/ 2 band with a high spectral resolution of 0.0044 \im 
(1 cm -1 ) with the intensity divided by 5, and at the resolution of 0.065 fim corresponding 
to the Spitzer IRS SL2 spectrometer. Ortho and para lines are indicated at the top, with 
the arrows showing the strongest lines for both spin species. Calculations pertain to a 0.075' 
field of view radius centered on the nucleus position, (^(H^O) = 5 x 10 28 molec. s _1 , = 
1.760 AU, A = 1.409 AU, v exp = 0.8 km s" 1 , T kin = 40 K, and OPR = 3. 
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Fig. 3. — Model fits to the spectra of C/2003 K4 (LINEAR) extracted along the slit positions 
(see § m Fig. [[]). Capital letters on the top-left corners correspond to the labels defined in 
the caption of Fig. [TJ Data are shown in black with pipeline derived errorbars excluding SL2 
fringe uncertainty (see text § 12.11) . with the model fits (in red) superimposed. Model fitting 
was performed with the rotational temperature T rot and the ortho-to-para ratio taken as free 
parameters. The derived water band intensities are given in Table [TJ Derived T rot (K) and 
OPRs are given in Tables [TJ and [2] (see § 13.21) . The underlying continuum, described by a 
polynomial of degree 5-6, was also fit simultaneously. 
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Fig. 4.— Model fits to the spectra of C/2003 K4 (LINEAR) (see §E1 Fig. QJ. Same as Fig. EJ 
with the fitted continuum background subtracted. Data are shown in black with errorbars, 
with the model fits (in red) superimposed. Errorbars include here SL2 fringe uncertainty 
(see § 12. ip . The residuals are shown on the bottom. Model fitting was performed with 
the rotational temperature T rot and the ortho-to-para ratio taken as free parameters. The 
derived values are given in Tables [1] and [2] (see § 13. 2j) . 



-27- 




Fig. 5. — Emergent model intensity of the water z/ 2 band as a function of rotational tem- 
perature, T rot . The synthetic water spectra were generated at the spectral resolution corre- 
sponding to the Spitzer IRS SL2 spectrometer (0.065 /xm). Calculations pertain to a 0.075' 
field of view radius centered on the nucleus position, Q(H 2 0) = 1 x 10 28 molec. s -1 , r h = 
1.760 AU, A = 1.409 AU, v exp = 0.8 km s _1 , with T rot equal to 20 K (solid black line), 30 K 
(dotted black line), 50 K (solid red line), 70 K (dotted red line), and 90 K (solid blue line). 
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Fig. 6. — Rotational temperature as a function of offset with respect to nucleus position in 
comet C/2003 K4 (LINEAR). Data points (black dots with errorbars) correspond to the slit 
extractions shown in Fig [TJ Curves show the rotation temperature extracted from synthetic 
v 2 spectra computed with a kinetic temperature T kin of 30, 40, 50, and 100 K (blue, black, red 
and green lines, respectively) and different models of the electron density (described by the 
x ne parameter): x ne = 1 (plain lines), x ne = 0.5 (dotted lines), x ne = 0.2 (dashed lines). For 
clarity, results obtained with x ne = 0.5 are only plotted for T kin = 30 K; for T kin = 100 K, 
we only show the model output corresponding to x ne = 0.2. 
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Fig. 7. — Emergent model intensity of the water v 2 band as a function of the ortho-to- 
para (OPR) ratio. OPRs of 2.0 (solid red line), 2.5 (dotted black line), and 3.0 (solid black 
line) are depicted. The synthetic water spectra were generated at the spectral resolution 
corresponding to the Spitzer IRS SL2 spectrometer (0.065 /xm). Calculations pertain to a 
0.075' field of view radius centered on the nucleus position, Q(B.20) — 1 x 10 28 molec. s _1 , 
r h = 1.760 AU, A = 1.409 AU, v exp = 0.8 km s~\ and T rot = 30 K. 
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Fig. 8.— Model fits to the spectrum of C/2003 K4 (LINEAR) at 7.2" offset from the 
nucleus. Data are shown in black dots with errorbars that include noise from spectral 
fringing (see § 12. II) . Synthetic spectra with OPR = 2.47 and OPR = 3 are shown in red and 
in black dotted line, respectively. The former fit was obtained with the OPR set as a free 
parameter, in contrast to the second fit. The retrieved T ro t are 30.7±3.2Kand 30.5±3.3K, 
respectively. The residual spectrum from the former fit is shown in red on the bottom, on 
which are su perimposed in arbitr ary units an interstellar PAH spectrum typical of class A 



sources fromlPeeters et al.l (120021 ) (blue spectrum), and a model of carbonate emission from 



Lisse et al.l (120061 ) (green spectrum) 
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Fig. 9. — Total band intensity between 5.8 and 7 fim in C/2003 K4 (LINEAR) as a function 
of offset with respect to nucleus position. Data points (black dots with errorbars) correspond 
to the slit extractions shown in Fig [TJ Squares show expected water z/ 2 band intensities for 
Q(H 2 0) = 2.43 x 10 29 molec. s"\ v exp = 0.8 km s" 1 and /3(H 2 0) = 1.6 x 1(T 5 s^ 1 at r h = 
1 AU. 



-32 - 



i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 



C/2003 K4 (LINEAR) _ 
[1.8AU] 




5.0 5.5 6.0 6.5 7.0 7.5 8.0 

Wavelength (am) 



Fig. 10. — A comparison of the v 2 water bands in Oort Cloud comets C/2003 K4 (LINEAR), 
C/2004 Bl (LINEAR), and Jupiter-family comets 71P/Clark and 9P/Tempel 1 detected with 
the Spitzer IRS. The data have not been scaled or offset in flux with respect to each other 
and the spectra do not have the continuum emission removed. Our synthetic model spectrum 
suggests that water dominates emission in-excess of the continuum at wavelengths between 
5.7 and 7 jum. 
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Fig. 11. — Background-subtracted 7-7.6 fim spectrum of C/2003 K4 (LINEAR) showing 
the detection of a weak emission feature at 7.3 fim. This spectrum has been obtained by 
averaging slit extractions 0-27, 0-22, 1-27, and 1-22. The uncertainty plotted are the error 
from extractions of the pipeline reduced individual BCDs and do not include any small excess 
noise arising from spectral fringing (see § 12. ip . 



